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RADAR-ESM  CORRELATION  DECISION  PROCEDURE 


INTRODUCTION 

The  problem  considered  in  this  report  is  how  to  correlate  an  Electronic  Support  Measures  (ESM) 
signal  with  one  or  none  of  m  possible  radar  tracks.  Previous  work  [1,2]  was  concerned  only  with  the 
angular  measurements  (i.e.,  it  excluded  all  classification  information  pertaining  to  the  type  of  radiating 
platform)  and  always  assumed  that  the  ESM  signal  went  with  one  of  two  radar  tracks.  Furthermore,  it 
was  assumed  that  the  number  of  measurements  associated  with  both  radar  tracks  was  the  same.  Refer¬ 
ence  1  was  concerned  with  several  least-squared  error  discriminants  and  compared  their  performance 
via  simulation.  Reference  2  was  concerned  with  the  problem  that  radar  and  ESM  measurements  are 
not  coincident  in  time,  so  that  measurements  from  one  or  both  sources  must  be  smoothed  to  common 
times.  The  present  report  again  is  only  concerned  with  angular  measurements,  but  it  considers  the 
additional  complexities  of  more  than  two  tracks  and  unequal  numbers  of  measurements.  The  possibil¬ 
ity  of  the  ESM  signal  not  belonging  to  any  one  of  the  radar  tracks  will  also  be  considered. 

LEAST-SQUARED  ERROR  DEVELOPMENT 

Assume  that  n  independent  and  unbiased  Gaussian-distributed  ESM  measurements  are  made  with 
a  variance  of  c r2.  Coleman  [2]  has  shown  that  the  probability  of  correct  association  is  only  slightly 
dependent  on  the  radar  measurement  accuracy  if  the  radar  accuracy  is  good  compared  to  the  ESM  accu¬ 
racy.  Thus,  for  simplicity  we  will  initially  assume  that  the  radar  measurements  have  zero  measurement 
error  (the  case  of  nonperfect  radar  measurements  is  considered  in  the  last  section),  and  we  will  further¬ 
more  assume  that  the  radar  measurements  are  coincident  in  time  with  the  ESM  measurements.  Then 
the  probability  of  error  can  be  minimized  [1]  by  choosing  the  radar  track  with  the  smallest  squared 
error, 

dj  =  £  [0,  -  0,  O’)]2.  7=1 . m.  (1) 

/-I 

where  9,  is  the  /  th  ESM  measurement  and  /3,  (J)  is  the  true  position  of  the  j  th  radar  target  at  the  time 
of  the  /  th  ESM  measurement.  This  decision  rule  assumes  that  the  m  possible  associations  are,  a  priori, 
equally  likely. 

Let  us  now  consider  the  problem  of  unequal  numbers  of  radar  measurements.  In  its  simplest 
form,  we  are  given  two  ESM  measurements,  9X  and  02>  two  positions  for  radar  track  1,  /?f  (1)  and 
P 2  0),  but  only  one  position  for  radar  track  2,  /3j  (2).  Since  the  corresponding  hypothesis- testing  prob¬ 
lem  is 

//,:  £(#,)  =  /J,  (1)  and  £|«2)  =  02<1) 

H2'.  E{9\\  =  /3|  (2)  and  £{02)  =  anything. 

no  uniform  most  powerful  test*  exists,  that  is,  no  optimal  test  exists  with  respect  to  the  probability  of 
correct  decision. 


Manuscript  submitted  on  January  c),  (MSI 

A  uniform  most  powerful  test  maximizes  the  probability  of  detection  for  any  point  m  the  composite  hypothesis  The  ah 
scncc  of  such  a  test  indicates  that  there  is  no  optimal  test  and  subopt imal  tests  must  be  used 
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One  way  of  proceeding  is  the  generalized  likelihood  ratio  test.  This  implies  substituting  maximum 
likelihood  estimates  for  the  unknown  parameters.  If  this  were  done,  any  unknown  fi,  (j)  would  simply 
be  replaced  by  0,.  Thus,  if  there  were  n}  radar  measurements  for  the  j  th  radar  track,  the  decision  rule 
(1)  would  reduce  to  choosing  the  smallest  squared  error 

dj  *  £  1*.  ~  &  (/)12.  j  *  1,  •  •  ■ .  w.  (2) 

/-I 

The  decision  rule  given  by  (2)  biases  the  decision  in  favor  of  the  radar  track  with  the  smallest  number 
of  measurements.  However,  (2)  is  not  the  only  option.  One  obvious  variation  of  it  is 

dj  -  —  ^  l 0/  -  Pj  0)]2.  j  —  1 . nt,  (3) 

ni  t-\ 

and  since  the  expected  value  of  the  summand  is  a2  when  j  is  the  correct  choice,  another  variation  is 

dj  (n  -  t\j)<j2  +  ^  [0,  -  Pi  (y)]2,  (4) 

i-i 

where  n  =  max  {nx,  n2 . nm }. 

To  consider  the  applicability  of  these  various  discriminants,  consider  the  following  two  extreme 
examples  in  which  there  are  100  samples  from  radar  track  1  and  only  one  sample  from  radar  track  2.  In 
the  first  example,  let  the  corresponding  squared  errors  [as  given  by  (2)1  be  90  a-2  and  0.8  or2,  respec¬ 
tively.  In  this  case  discriminants  (2)  and  (3)  choose  track  2  and  discriminant  (4)  chooses  track  1 . 
Given  the  fact  that  10(kr2  is  the  expected  value  of  the  squared  error  (if  the  correct  association  is  with 
radar  track  1),  we  believe  the  decision  should  be  made  to  associate  the  ESM  signal  with  track  1.  This  is 
because  the  probability  of  obtaining  a  squared  error  of  at  least  90cr2  from  100  measurements  from  the 
correct  track  is  near  1,  and  this  is  larger  than  the  probability  of  obtaining  a  squared  error  of  at  least 
O.8o*2  from  one  measurement  from  the  correct  track.  In  the  second  example,  the  corresponding 
squared  errors  are  104o-2  and  3<r2.  In  this  case  discriminants  (2)  and  (4)  choose  track  2  and  discrim¬ 
inant  (3)  chooses  track  1.  The  same  probabilistic  argument  again  indicates  that  track  1  is  the  correct 
track.  Thus,  in  the  two  examples  none  of  the  three  discriminants  consistently  yielded  the  "correct11 
track.  Rather,  we  believe  that  the  decision  should  be  based  on  the  largest  probability  of  obtaining  a 
squared  error  greater  than  the  measured  value,  given  that  the  radar  target  under  consideration  is  the 
correct  one. 

PROBABILITY  STATISTICS 

The  suggested  procedure  is  as  follows:  for  each  radar  track,  calculate  the  probability  of  obtaining 
at  least  the  observed  squared  error  under  the  assumption  that  the  ESM  signal  comes  from  the  radar 
track.  This  is  accomplished  by  first  normalizing  the  squared  error;  i.e., 

s,  =  £  b,  -  0,  0)12/<r2.  (5) 

I- I 

Then,  if  the  association  is  the  correct  association  (i.e.,  E{0j]  =  P,(j)y  i  *  I . n} ),  Sj  has  a  chi- 

square  density  with  rtj  degrees  of  freedom.  Thus,  the  desired  probability  P}  is  given  by 

-  Pj  =  Pr  [z  >  S;).  where  z  —  x2 

It  is  sho^n  in  the  appendix  that 

Prizes}-/*,  (6) 

where  /*  is  obtained  by  iteration, 

/*  “  fk  +  /*-!.  (7) 

F*-  Fk^S/2k. 
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where  the  initial  values  are 

#0-  fo«  *-s72  (8) 

and  the  even  integer  n  is  related  to  k  by 


Now  define  Pm ax  as  the  largest  P,  and  Pnext  as  the  second  largest  P}  in  the  set  of  P/s,  j  =  1  ...  m. 
From  this  set  of  P/s,  one  of  five  possible  decisions  is  made: 

(1)  ESM  signal  goes  with  track  corresponding  to  Pmax; 

(2)  ESM  signal  probably  goes  with  track  corresponding  to  Pmax; 

(3)  ESM  signal  probably  goes  with  some  track; 

(4)  ESM  signal  probably  does  not  go  with  any  track;  or 

(5)  ESM  signal  does  not  go  with  any  track. 

We  arrive  at  these  decisions  by  defining  three  probability  thresholds,  TH  (high),  TSi  (middle),  and  TL 
(low),  and  a  probability  ratio  R.  The  corresponding  decision  rules  are 

Category  1 .  Pmax  ^  atld  Pmax  ^  R  Pnexl i 

Category  2.  7//  >  Pmax  ^  and  Pmax  ^  ^  Pnext* 

Category  3:  P,™  >  TM  but  Pmax  <  P  Pnexf  \ 

Category  4:  TM  >  Pmax  >  TL  \  and 

Category  5:  TL  >  P^x 

Furthermore,  if  for  any  P,  both  P,  <  TL  and  R  P,  <  Pmax  are  true,  the  /  th  radar  track  is  dropped  from 
further  consideration  for  association  with  the  ESM  signal.  It  should  be  noted  that  if  the  number  of 
radar  measurements  is  the  same  for  each  track,  the  outlined  procedure  reduces  to  picking  the  radar 
track  with  the  minimum  squared  error.  The  problem  which  will  now  be  considered  is  how  to  set  the 
three  thresholds,  Tt ,  TMy  and  TH. 

THRESHOLD  SETTINGS 

The  easiest  threshold  to  set  is  TL.  This  threshold  determines  the  probability  that  the  correct  radar 
track  (i.e.,  the  one  associated  with  the  ESM  signal)  will  be  incorrectly  rejected  from  further  considera- 
tion.  The  key  element  in  setting  the  threshold  is  noting  that  the  probability  P,  for  the  correctly  associ¬ 
ated  radar  track  is  uniformly  distributed  between  zero  and  one.  This  can  be  shown  as  follows:  First, 
note  that  Pj  is  given  by 

Pj  =  F  ($)  &  f  ~  pxAz)  dz.  (9) 

where  p  2  (*)  is  the  probability  density  function  of  a  chi-square  random  variable  with  n,  degrees  of  free¬ 
dom.  The  value  of  P}  is  obviously  bounded  between  zero  and  one.  The  cumulative  distribution  of  P, 
can  then  be  written  as 

Pr  {Pj  ^  Pi  -  Pr\F(S,)  ^  P|~  Pr [S,  >  F  1  (P)}.  (10) 
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where  the  direction  of  the  inequality  changes  because  F  (•)  is  a  monotonic  decreasing  function  of  its 
argument.  When  target  j  is  the  correct  association,  Sj  [defined  by  (5)]  has  the  density  p^ASj)  men¬ 
tioned  in  Eq.  9.  Therefore, 

Pr  {Sj>z)  =  F(z)  (ID 

from  (9).  Substituting  (11)  into  (10)  yields 

Pr  [Pj  <  P)  =  F[r'(P))  =  P. 

and  thus,  Pj  is  uniformly  distributed.  Consequently,  if  one  desires  a  rejection  rate  of  PR ,  one  can 
obtain  this  by  setting  TL  =  PR. 

In  a  similar  vein  of  thought  the  threshold  TH  is  set  by  PfA ,  the  probability  of  falsely  associating  a 
track  with  an  ESM  signal  when  the  signal  does  not  belong  with  any  of  the  tracks.  Then  Sj  [defined  by 
(5)1  has  a  noncentral  chi-square  density  with  noncentrality  parameter 

Fk  =  £  Iflflil  -  P,  (j))1/*2.  (12) 

J  1-1 

Then  corresponding  to  (11)  one  has 

Pr  {Sj  >  2}  =  Fk  (2),  (13) 

where  Fk  (■)  is  the  noncentral  chi-square  distribution.  Substituting  (13)  into  (10)  we  obtain 

Pr  [Pj  <  P]  -  FK  [r'(P)  1.  (14) 

The  threshold  TH  is  now  determined  by 

PFa  =  Pr  [Pj  >  TH)=  1  -  Fklr'(TN)l  (15) 

Unfortunately,  (15)  cannot  be  inverted  to  yield  an  analytic  expression  for  TH. 

Obviously,  the  threshold  TH  is  a  function  of  which  is  related  to  the  difference  between  the 
true  (ESM)  position  and  the  radar  track  under  consideration.  This  difference  was  denoted  by  /lc. 

m  =  £{e,)  -j3,0). 

and  the  threshold  TH  was  found  for  arbitrarily  chosen  values  of  \i  by  simulation  techniques.  The 
results  for  PhA  =  0.01  and  fi=\.0a  and  1.5<r  are  shown  in  Fig.  1.  For  each  value  of 

nj  (2,4,6 .  100)  3000  repetitions  were  performed.  The  threshold  is  near  one  for  small  values  of  n 

and  drops  to  zero  as  n  approaches  infinity.  Since  the  threshold  is  still  high  for  large  n's  when  fi  =  l  .Ocr, 
the  curve  for  fx  =  1.5a*  was  used  in  later  simulations.  Of  course,  the  value  that  should  be  used  in  a  real 
system  would  depend  on  the  expected  target  separation  divided  by  the  measurement  accuracy  in  a  real¬ 
istic  scenario. 

The  middle  threshold  divides  the  "tentative"  region  into  a  tentative  track  region  and  a  tentative  no 
track  region.  The  rationale  in  setting  the  threshold  is  to  set  both  of  these  probabilities  equal  for  a  par¬ 
ticular  separation;  i.e., 

Pr  [Pj  <  Tm  I  correct  match)  -  Pr{P,  ^  TM  \  incorrect  match). 

Thus  the  equation  specifying  is 

TM=  1  -FKj[rlw„)\. 

Again,  the  threshold  is  a  function  of  and  must  be  found  by  simulation  techniques.  The  previous 
simulation  used  to  find  Th  was  also  used  to  find  Ty ;  the  results  are  shown  in  Fig.  2.  Again,  the  curve 
for/*  -  l.Sfr  was  used  for  further  simulation  results. 
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Tig  2  —  Middle  threshold  vs  number 
of  samples  for  two  different  target 
separations 


The  only  parameter  that  remains  to  be  set  is  the  probability  ratio  R,  which  is  used  to  inhibit  a  firm 
track  decision  if  is  not  greater  than  R  times  No  methodology  has  been  developed  for  setting 
R.  In  the  simulation  R  has  been  set  to  9.  (Some  of  the  simulation  results  suggest  that  this  value  may 
be  too  large.) 

SIMULATION  RESULTS 

All  simulations  involved  an  attempt  to  associate  an  ESM  signal  with  one  of  three  possible  radar 
tracks  whose  angular  positions  are  known  to  be  #*t,  *i2,  and  that  is,  0,  (./)  »  M)  for  all  measure¬ 
ments.  First,  two  ESM  angles  of  arrival  were  generated  from  a  zero-mean  Gaussian  density  with  unit 
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variance.  The  discriminant  probabilities  P),  P;.,  and  P*  were  then  calculated  tor  n  «  2  using  the  itera¬ 
tive  relation  given  by  (6),  (7),  and  (8).  The  highest  probability  Pmax  was  found  and  compared  to  the 
various  thresholds  to  determine  to  which  of  the  five  possible  categories  the  ESM  signal  should  be 
assigned.  In  all  simulations  the  lower  threshold  was  set  to  0.00 1,  the  upper  threshold  was  set  by  the 
M  =  1.5<r  curve  in  Fig.  1,  the  middle  threshold  was  set  by  Uicm  =  l  5<r  curve  in  Fig.  2,  and  the  proba¬ 
bility  ratio  was  set  at  9.0.  Next,  two  more  HSM  measurements  were  generated  and  the  procedure  just 
described  was  repeated  for  n  —  4  measurements.  This  addition  of  two  ESM  measurements  with  each 
repetition  continued  until  50  ESM  measurements  were  generated.  The  whole  set  of  measurements 
comprises  a  single  radar-ESM  correlation  trial.  To  obtain  significant  statistics,  1000  repetitions  were 
generated  for  each  case. 

For  case  1,  mi  »  0.  M2  “  1,  and  M  3  *  -  1;  the  results  of  1000  repetitions  are  shown  in  Table  1. 
Initially,  when  only  2  ESM  measurements  are  present,  approximately  80%  of  the  ESM  signals  are  tenta¬ 
tively  associated  with  no  radar  track.  The  number  of  tentative  associations  with  no  track  falls  below 
10%  with  16  ESM  measurements,  and  it  eventually  falls  below  2%  with  50  ESM  measurements.  The 
number  of  cases  in  which  a  firm  "no  track"  decision  is  made  is  always  less  than  0.5%.  While  originally 
80%  of  the  ESM  signals  are  tentatively  associated  with  unknown  tracks,  it  is  not  until  one  receives  26 
ESM  measurements  that  50%  of  the  ESM  signals  are  firmly  correlated  with  the  correct  track,  track  1. 
This  is  caused  by  two  facts:  First,  the  high  threshold  is  0.29  for  26  ESM  measurements.  Second,  the 
difference  between  71%  (which  is  the  percentage  above  the  high  threshold)  and  the  actually  obtained 
50%  firm  decisions  is  due  to  the  probability  ratio,  R  =  9. 


Table  1  —Simulation  Results  for  Five-Decision  Association 
Rule  as  a  Function  of  the  Number  of  ESM  Measurements  (Target 
separations  areMi  0.  M2  *r  1.  and  M3  s 


Number  of  HSM 

Firm  Correlation 

Tentative  Correl 

ation 

Tentative 

Unknown  Track 

Tentative 

No  Track 

Firm 

No  Track 

Measurements 

Track  1 

Track  2 

Track  3 

Track  1 

Track  2 

Track  3 

2 

0 

0 

0 

0 

3 

10 

818 

169 

0 

4 

0 

0 

0 

0 

6 

6 

786 

201 

1 

6 

0 

0 

0 

0 

2 

4 

804 

190 

0 

8 

0 

0 

0 

0 

2 

4 

840 

154 

0 

10 

0 

0 

0 

23 

2 

2 

836 

137 

0 

12 

0 

0 

0 

106 

l 

I 

787 

105 

0 

14 

23 

0 

0 

188 

2 

0 

674 

113 

0 

16 

(00 

1 

0 

240 

0 

0 

562 

96 

1 

18 

191 

l 

0 

232 

0 

1 

484 

91 

0 

20 

279 

0 

1 

218 

0 

0 

421 

81 

0 

22 

379 

1 

0 

192 

l 

0 

355 

71 

1 

24 

433 

0 

0 

191 

l 

0 

312 

62 

1 

26 

512 

0 

0 

166 

1 

0 

264 

54 

3 

28 

597  1 

0  1 

0 

149 

0 

0 

206 

45 

3 

30 

660 

0 

0 

123  ! 

o 

0 

172 

42 

3 

32 

705 

0 

0 

112 

0 

0 

141 

39 

3 

34 

765 

0 

0 

88 

0 

0 

106 

38 

3 

36 

796 

0 

0 

70 

0 

0 

95 

36 

3 

38 

828 

0 

0 

55 

0 

0 

84 

29 

4 

40 

858 

0 

0 

40 

0 

0 

69 

29 

4 

42 

874 

0 

0 

41 

0 

0 

58 

22 

5 

44 

895 

0 

0 

38 

0 

0 

40 

23 

4 

46 

912 

0 

0 

27 

0 

0 

35 

21 

5 

48 

933 

0 

0 

18 

0 

0 

24 

21 

4 

50 

936 

0 

0 

20 

0 

0 

22 

18 

4 
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For  case  2,  with  fi\  *  2.  ~  K  and  /la3  =  -  1,  the  results  of  1000  repetitions  are  shown  in  Table 

2.  In  this  case  the  ESM  signal  does  not  belong  to  any  radar  track  but  lies  between  radar  tracks  2  and  3. 
Initially,  about  75%  of  the  ESM  signals  are  tentatively  associated  with  radar  tracks,  and  after  8  ESM 
measurements  this  number  is  reduced  to  50%.  It  is  not  until  one  has  slightly  less  than  50  ESM  meas¬ 
urements  that  90%  of  the  ESM  signals  are  either  tentatively  or  firmly  declared  not  to  be  associated  with 
any  radar  track.  It  should  be  noted  that  while  initially  75%  of  the  ESM  signals  were  tentatively  associ¬ 
ated  with  radar  tracks,  at  no  time  was  the  number  of  firm  associations  more  than  5%. 


Table  2— Simulation  Results  for  Five-Decision  Association 
Rule  as  a  Function  of  the  Number  of  ESM  Measurements  (Target  separations 
are  ~  2,  —  I,  and  ixy  —  —  1) 


Number  of  ESM 

Measurements 

Firm  Correlation 

Tentative  Correlation 

Tentative 

Unknown  Track 

Tentative 

No  Track 

Firm 

No  Track 

Track  1 

Track  2 

Track  3 

Track  1 

Track  2 

Track  3 

2 

0 

o 

U 

0 

0 

183 

545  ^ 

261 

0 

4 

0 

0 

11 

0 

65 

156 

388 

378 

2 

6 

0 

n 

16 

0 

141 

158 

206 

468 

0 

8 

0 

n 

24 

0 

152 

143 

169 

500 

1 

10 

0 

20 

11 

0 

130 

153 

124 

559 

3 

12 

0 

19 

14 

0 

144 

166 

109 

539 

9 

14 

0 

25 

18 

0 

138 

150 

7! 

579 

19 

16 

0 

20 

19 

0 

142 

127 

58 

601 

33 

18 

0 

24 

16 

0 

133 

113 

45 

627 

42 

20 

0 

25 

14 

0 

115 

no 

40 

635 

61 

22 

0 

28 

17 

0 

102 

107 

29 

628 

89 

24 

0 

23 

18 

0 

112 

102 

23 

605 

117 

26 

0 

26 

23 

0 

93 

89 

18 

608 

143 

28 

0 

24 

21 

0 

87 

76 

11 

603 

178 

30 

0 

25 

24 

0 

78 

66  1 

9 

582 

216 

32 

0 

24 

23 

0 

62 

53 

8 

599 

231 

34 

0  ! 

24  | 

20 

0 

60 

54 

H 

555 

276 

36  1 

0  1 

24 

21 

0 

54 

50 

6 

535 

310 

38 

0 

25 

21 

0 

52 

45 

6 

502 

349 

40  | 

0 

30 

20 

0 

47 

40 

5 

462 

396 

42  | 

0 

29 

21 

0 

40 

32 

2 

428 

448 

44 

0 

27 

21 

0 

32 

27 

4 

393 

4% 

46 

0 

28 

21 

0 

31 

20 

2 

375 

523 

48 

0 

29 

22 

0 

23 

20 

2 

343 

561 

50 

0 

31  1 

21 

0 

24 

20 

I 

304 

1  599 

In  both  cases,  the  results  improve  dramatically  as  the  number  of  ESM  measurements  increases. 
However,  the  probability  of  correct  decision  does  not  reach  90%  until  about  50  measurements  are 
received.  It  should  be  noted,  however,  that  the  situation  encountered  in  cases  1  and  2  is  rather 
difficult:  there  is  an  incorrect  track  one  standard  deviation  away  from  the  true  track.  Consequently, 
cases  3  and  4  were  generated.  In  case  3,  mi  =  0,  fi2  “  2,  and  /iy  =*  -  2,  and  the  results  are  given  in 

Table  3.  If  we  compare  Tables  1  and  3,  while  the  number  of  tentative  no  track  decisions  is  slightly 

smaller  for  Table  3  (2 o*  separation),  the  major  difference  is  that  the  ESM  signals  are  associated  with 
track  1  very  quickly.  In  particular,  the  number  of  associations  with  "unknown"  tracks  falls  to  0.2%  after 

only  10  ESM  measurements.  In  case  4,  *  4,  “  2,  and  M3  “  ~  2,  and  the  results  are  given  in 

Table  4.  Here,  after  10  ESM  measurements,  90%  of  the  ESM  signals  are  correctly  designated  as  not 
associating  with  any  radar  track,  and  at  22  ESM  measurements,  all  ESM  signals  are  correctly  designated 
as  not  associating  with  any  radar  track.  Thus,  when  the  nearest  incorrect  radar  track  is  two  standard 
deviations  away  from  the  correct  track,  correct  decisions  are  arrived  at  very  rapidly. 

One  of  the  mayor  reasons  for  basing  the  decision  on  the  probability  statistic  P,  is  the  resulting  abil¬ 
ity  to  handle  a  different  number  of  radar  measurements  for  each  radar  track.  To  test  this  behavior, 
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Table  3— Simulation  Results  for  Five-Decision  Association 
Rule  as  a  Function  of  the  Number  of  HSM  Measurements 
(Target  separations  are/X|  =  0,  =  2,  and/13  =  ~2) 


Number  of  FSM 

Measurements 

Firm  Correlation 

Tenatative  Correlation 

Tentative 

Unknown  Track 

Tentative 

No  Track 

Firm 

Track  1 

Track  2 

Track  3 

Track  1 

Track  2 

Track  3 

No  Track 

2 

“26" 

0 

0 

371 

10 

14 

339 

240 

0 

4 

83 

0 

1 

536 

4 

4 

99 

272 

1 

6 

139 

0 

0 

576 

0 

2 

34 

249 

0 

8 

216 

0 

0 

586 

0 

0 

12 

184 

2 

10 

276 

0 

0 

565 

0 

0 

2 

156 

1 

12 

334 

0 

0 

539 

0 

0 

1 

125 

I 

14 

409 

0 

0 

469 

0 

0 

0 

120 

2 

16 

458 

0 

0 

439 

0 

0 

0 

100 

3 

18 

520 

0 

0 

384 

0 

0 

0 

95 

1 

20 

581 

0 

0 

335 

0 

0 

0 

82 

2 

22 

646 

0 

0 

278 

0 

0 

0 

75 

l 

24 

683 

0 

0 

251 

0 

0 

0 

64 

2 

26 

729 

0 

0 

211 

0 

0 

0 

59 

1 

28  1 

777  1 

0 

0  j 

173 

0 

0 

0 

49 

1 

30 

813 

0 

0  1 

139 

0 

0 

0 

47 

1 

32 

829 

0  ! 

0 

129 

0 

0 

0 

41 

1 

34 

860 

0 

0 

100 

0  1 

0 

0 

39 

1 

36  1 

885 

0  I 

0  1 

77 

0 

0 

0 

37 

1 

38  ! 

907 

0 

0 

61 

0 

0 

0 

30 

2 

40 

925 

0 

0 

43 

0 

0 

0 

30 

2 

42 

931 

0 

0 

43 

0 

0 

0 

23 

3 

44  | 

936 

0 

0 

38 

0 

0 

0 

24 

2 

46  1 

948 

0 

0 

27 

0 

0 

0 

22 

3 

48 

958 

0 

0 

18 

0 

0 

0 

22 

2 

50 

959 

0 

0 

20 

0 

0 

0 

19 

2 

Table  4— Simulation  Results  for  Five-Decision  Association  Rule  as  a  Function  of  the  Number  of 
ESM  Measurements  (Target  separations  are  /lx  t  =  4.  =  2,  and  =  -2) 
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cases  5  and  6  were  generated.  Both  cases  are  similar  to  case  3  in  that  Mi  =*  0,  /a  2  =  2,  and  fx}=  -  2. 
The  difference  is  that  in  case  5  one  starts  off  with  12  ESM  measurements  from  radar  track  1  but  with 
only  2  radar  measurements  each  from  radar  tracks  2  and  3,  and  in  case  6  one  starts  off  with  12  ESM 
measurements  and  12  radar  measurements  each  from  radar  tracks  2  and  3,  but  only  2  radar  measure¬ 
ments  from  radar  track  1.  The  results  for  cases  5  and  6  are  given  in  Tables  5  and  6,  respectively.  If  we 
compare  Tables  3,  5,  and  6,  the  following  conclusions  can  be  drawn: 

(1)  The  number  of  associations  with  unknown  tracks  is  primarily  determined  by  the  number  of 
measurements  on  -incorrect  tracks;  i.e.,  after  a  reasonable  number  of  measurements  is  obtained  on 
incorrect  tracks,  they  can  be  recognized  as  not  being  associated  with  the  ESM  signal.  (Note:  12  ESM 
measurements  in  Table  3,  22  ESM  measurements  in  Table  5,  and  12  ESM  measurements  in  Table  6  all 
correspond  to  12  radar  measurements  on  tracks  2  and  3.) 

(2)  The  number  of  tentative  "no  track"  decisions  is  mainly  a  function  of  the  number  of  measure¬ 
ments  on  the  correct  track,  track  1.  (The  12,  12,  and  22  ESM  measurements  in  Tables  3,  5,  and  6, 
respectively,  all  correspond  to  12  measurements  on  track  1.) 

(3)  The  numbers  of  firm  and  tentative  decisions  to  associate  the  ESM  signal  with  the  correct 
track,  track  1,  are  again  mainly  functions  of  the  number  of  radar  measurements  on  the  correct  track. 

Thus,  it  appears  that  the  proposed  method  yields  good  results  when  the  number  of  measurements  asso¬ 
ciated  with  each  radar  track  is  different. 


Table  5— Simulation  Results  for  Five-Decision  Association  Ride  as  a  Function 
of  the  Number  of  ESM  Measurements  (Target  separations  are  mi  *  0.  m:  “  2,  and  M3  -  -2; 
initial  number  of  measurements  on  each  radar  track  is  12,  2,  and  2,  respectively) 


Number  of  ESM 

Measurements 

Firm  Correlation 

Tentative  Correlation 

Tentative 

Unknown  Track 

Tentative 

No  Track 

Firm 

No  Track 

Track  1 

Track  2 

Track  3 

Track  1 

Track  2 

Track  3 

12 

225 

0 

0 

169 

8 

6 

455 

137 

*0  1 

14 

365 

0 

0 

375 

1 

0 

140 

119 

0 

16 

455 

0 

0 

384 

0 

0 

52 

109 

0 

18 

496 

0 

0 

393 

0 

0 

12 

96 

3 

20 

588 

0 

0 

316 

0 

0 

3 

91 

2 

22 

629 

0 

0 

290 

0 

0 

1 

76 

4 

24 

671 

0 

0 

259 

0 

0 

1 

65 

4 

26 

728 

0 

0 

210 

0 

0 

1 

58 

3 

28 

779 

0 

0 

162 

0 

0 

0 

56 

3 

30 

802 

0 

0 

153 

0 

0 

0 

42 

3 

32 

838 

0 

0 

119 

0 

0 

0 

39 

4 

34 

858 

0 

0 

107 

0 

0 

0 

33 

2 

36 

889 

0 

0 

74 

0 

0 

0 

35 

2 

38 

902 

0 

0 

70 

0 

0 

0 

27 

1 

40  1 

923 

0 

0 

49 

0 

0 

0 

27 

1 

42 

929 

0 

0 

42 

0 

0 

0 

26 

3 

44 

940 

0 

0 

35  ! 

0 

0 

0 

22 

3 

46 

952 

0 

0  ! 

27 

0 

0 

0 

18 

3 

48  1 

964 

0  1 

0  I 

18  | 

0 

0 

0 

15 

3 

50 

968 

0 

o 

17 

0 

0 

0 

13 

2 

NONPERFECT  RADAR  MEASUREMENTS 

In  the  previous  development  it  was  assumed  that  the  radar  measurements  were  perfect.  We  will 
now  assume  that  the  radar  makes  noisy  measurements,  that  0,  (/)  is  the  predicted  azimuth  of  the  j  th 
radar  track  at  the  time  of  the  i  th  ESM  measurement  and  is  Gaussian  distributed,  and  that  the  means 
and  variances  of  0,  and  p(  (j)  are 
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Table  6— Simulation  Results  for  Five-Decision  Association  Rule  as  a  Function 
of  the  Number  of  ESM  Measurements  (Target  separations  are  (jl\  =  0,  =  2,  and  =  -2; 

initial  number  of  measurements  on  each  radar  track  is  2,  12,  and  12,  respectively) 


Number  of  ESM 

Firm  Correlation 

Tentative  Correlation 

Tentative 

Tentative 

Firm 

Measurements 

Track  1 

Track  2 

Track  3 

Track  1 

Track  2 

Track  3 

Unknown  Track 

No  Track 

No  T  rack 

12 

27 

0 

0 

650 

0 

r  o 

0 

321 

2 

14 

75 

0 

0 

658 

0 

0 

0 

266 

1 

16 

160 

0 

0 

610 

0 

0 

0 

230 

» 

18 

201 

0 

0 

601 

0 

0 

') 

198 

o  ; 

20 

292 

0 

0 

535 

0 

0 

0 

172 

I  ; 

22 

349 

0 

0 

507 

0 

0 

0 

140 

4 

24 

415 

0 

0 

451 

0 

0 

0 

131 

3 

26 

482 

0 

0 

393 

0 

0 

0 

123 

2 

28 

536 

0 

0 

349 

0 

0 

0 

113 

2 

30 

612 

0 

0 

288 

0 

0 

0 

98 

2  ! 

32 

660 

0 

0 

260 

0 

0 

0 

78 

!  2  : 

34 

672  ! 

0  ! 

o 

253 

o  ! 

o  1 

0 

73  1 

36 

727 

0 

0 

212 

0  , 

0  , 

0 

60 

1  i  i 

38  | 

758 

0  I 

0 

183  | 

0 

0 

0 

57  i 

? 

40 

804 

0 

0 

151 

0 

0 

0 

43  | 

2 

42  1 

828 

0 

0 

123 

0 

0 

0 

47 

i  2 

44 

862 

0 

0 

99 

0 

0 

0 

37  1 

2 

46 

883  j 

0 

0 

84 

0 

0 

0 

30 

3 

48 

906 

0  1 

0 

66 

0 

0 

0 

25  ! 

50  | 

916  | 

0  1 

0 

62  ! 

0  I 

0 

0 

20  | 

1 _ 2 _  I 

E  {0,}  =  E  [fij  (j)},  j  =  true  track, 
Var{0,)  =*  <T£t 
and 


Var{0,  (/))  = 

Now,  if  one  sets 

cr 2  =  c t\  -F  <r}r 

the  statistic  S}  given  by  (5)  still  has  a  chi-square  density  for  the  correct  association.  Thus,  it  follows 
that  all  the  thresholds  should  have  the  same  value,  whether  or  not  the  radar  measurements  are  nois>. 
There  will  be  some  difference  in  the  multitarget  performance  because  of  the  different  dependence 

between  the  squared  errors  {S,  :  j  —  1 . m\  due  to  the  radar  variances.  As  shown  in  [2],  the 

degree  of  smoothing  of  the  radar  data  (or  degree  of  dependency  between  radar  measurements)  can  also 
affect  performance.  However,  in  this  simulation  the  measurements  were  unsmoothed  (independent). 

To  show  that  the  performance  difference  is  small,  case  3  with  fi\  ~  0,  ^  :  2,  and  ft  j  *  -  2  was 

repeated  with  cr l  —  0.06  and  o*,2  =  0.94,  and  the  results  are  shown  in  Table  7.  Comparing  the  results 
in  Table  7  with  those  in  Table  3,  one  sees  very  little  difference.  Thus,  the  total  variance  is  the  impor¬ 
tant  quantity:  it  makes  little  difference  with  this  algorithm  how  the  uncertainty  is  divided  between  the 
radar  and  the  ESM  measurements. 

SUMMARY 

When  the  numbers  of  measurements  associated  with  different  radar  tracks  are  not  the  same,  the 
squared  error  association  criterion  should  not  be  used.  Rather,  the  decision  should  be  based  on  the 
conditional  probability  P,  of  obtaining  at  least  the  observed  squared  error  given  that  the  association  is 
the  correct  one.  Using  this  probability,  five  possible  decisions  were  formulated: 
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Table  7  — Simulation  Results  for  Five-Decision  Association  Rule  as  a  Function 
of  the  Number  of  KSM  Measurements  (Target  separations  are  /a  ,  -  0 .  m  ?  ”  2,  and  ^  ~2y 

measurement  accuracies  are  <r  /  =  0.06  and  cr  ^  =  0.94 


Number  of  FSM 

Firm  Correlation 

Tentative  Correlation 

I  email  ve 

TcntaiivtF 

:  Amp 

Measurements 

Track  1 

T  rack  2 

Track  3 

Track  1 

T  rack  2 

T  rack  3 

Unknown  1  rack 

No  1  rack^ 

2 

0 

l 

336 

3 

2  ' 

330 

303 

0 

4 

62 

0 

0 

549 

0 

2 

112 

275 

0 

6 

136 

0 

0 

583 

0 

0 

28 

252 

1 

8 

201 

0 

0 

566 

0 

0 

11 

221 

1 

to 

251 

0 

0 

559 

0 

0 

S 

185 

0 

12 

301 

0 

0 

535 

0 

0 

1 

163 

0 

14 

373 

0 

0 

477 

0 

0 

1 

149 

0 

16 

444 

0 

0 

418 

0 

0 

0 

138 

0 

18 

491 

0 

0 

401 

0 

0 

1 

106 

l 

20 

551 

0 

0 

346 

0 

0 

0 

103 

0 

22 

601 

0 

0 

315 

0 

0 

0 

84 

0 

24 

646 

0 

0 

292 

0 

0 

0 

62 

0 

2n 

695 

0 

0 

239 

0 

0 

0 

66 

0 

28 

753 

0 

0 

186 

0 

0 

0 

61 

i  o 

30 

799 

0 

0 

144 

0 

0 

0 

57 

0 

32 

831 

0 

0 

113 

0 

0 

0 

56 

1  0 

34  1 

863 

0 

0 

89 

0 

0 

0 

48 

0 

36  1 

872 

0  I 

0 

90 

0 

0 

0 

38 

0  1 

38  ! 

9(H) 

0 

0 

63 

0 

0 

0 

.P 

'  l) 

40 

922 

0 

0 

48 

0 

0 

0 

29 

;  1 

42 

937 

0 

0 

30 

0 

0 

0 

33  ; 

0 

44 

941 

0 

0 

30 

0 

0 

0 

29  | 

0 

46 

942 

0 

0  1 

31 

t) 

0  1 

0 

17  ! 

0  " 

48 

950 

0 

0 

23 

0 

0 

0 

2o  : 

1 

50 

059 

0 

0 

21 

0 

o  1 

0 

19  1 

1  1 

FSM  signal  goes  with  a  particular  track; 

FSM  signal  probably  goes  with  a  particular  track; 

FSM  signal  probably  goes  with  some  track; 

FSM  signal  probably  does  not  go  with  any  track;  and 
FSM  signal  does  not  go  with  any  track. 

These  decisions  are  reached  by  comparing  the  largest  Pt  to  various  thresholds  and  to  the  next  largest  l\. 
Various  simulations  were  run  in  an  attempt  to  quantify  the  performance  in  several  situations.  The  per¬ 
formance  was  generally  good  when  targets  were  separated  by  twice  the  FSM  measurement  accurac> 
(standard  deviation). 

There  are  several  problems  still  remaining  concerning  the  radar-FSM  correlation  algorithms. 
Among  these  are: 

U)  Should  one  continually  accumulate  the  squared  error  of  all  the  samples  or  should  one  weight1 
the  past  squared  errors  differently  after  a  hxed  number  of  samples  has  been  accumulated? 

(2)  If  one  has  declared  that  an  FSM  signal  belongs  to  a  particular  radar  track,  how  long  should* 
one  keep  this  association  if  a  new  radar  track  with  a  small  Sf  and  corresponding  large  P,  appears  ' 
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(3)  How  should  the  algorithms  be  modified  to  take  into  account  equipment  deficiencies  like  bias 
and  quantized  measurements? 

(4)  How  should  nonangular  information  ( i.e. ,  radar-track  parameters  such  as  velocity  and  eleva¬ 
tion  and  ESM  parameters  such  as  signal  type  and  equipment  identification)  be  included  in  the  discrim¬ 
inant  function? 

Future  work  will  be  concerned  with  answering  these  questions  and  testing  the  algorithms  with 
recorded  experimental  data. 
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Appendix 

CALCULATION  OF  CHI-SQUARE  PROBABILITY 

The  probability  that  z  is  greater  than  s  when  z  is  chi-squared  distributed  with  n  degrees  of  freedom 
is  given  by 


PrU  >  s]  =  js 


zinl 2)~\  e~zt 2 


-  Js  2^r(n/2)^ 

Assuming  that  n  is  an  even  integer  and  defining  A:  by  k  —  (n/2)  —  1,  we  get 

.  %  C°°  ike~zjl 

I,  -  Friz  »  >|  -  -^7*. 


Integrating  by  parts  yields 


which  can  be  simplified  to 


—zk  e~z ^ 


2*0c  -  1) !  ’ 


p-s/2 

L  = 

*  2kk'. 


Thus,  4  can  be  calculated  interatively  by 


4  **  Fk  +  4-i. 

*  -  TZ  Fk-u 


where 


/o- Fo~  e~s". 
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